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Local  Geoid  and  Gravity  Anomaly 
Predictions  Using  Point  Masses 


1.  INTRODUCTION 

In  order  to  represent  the  geoid  to  a  desired  resolution  using  spherical  har¬ 
monics.  the  necessary  number  of  potential  coefficients  is  exceedingly  large.  For 
example,  the  resolution  corresponding  roughly  to  1° (about  111  km)  requires  a  set 
of  coefficients  complete  through  degree  and  order  (180,  180)  olving  almost 
33,  000  coefficients.  To  produce  a  resolution  twice  as  fine,  that  is,  about  0.  5°, 
four  times  as  many  coefficients  are  required.  If  enormous  quantities  of  data  were 
available  all  over  the  globe  so  that  such  sets  of  coefficients  would  be  theoretically 
solvable  in  a  least  squares  adjustment,  the  matrix  of  normal  equations  in  the  1°  case 
would  have  dimensions  of  almost  33,  000  X  33,  000.  The  difficulties  connected  with 
the  actual  solution  of  this  system  would  be  formidable.  Thus,  even  to  describe  the 
geoid  over  a  relatively  small  region,  such  as  the  North  Atlantic,  would  be  extremely 
wasteful  if  not  impossible. 

The  use  of  local  parameters  called  "point  masses"  has  been  envisioned  as  a 
technique  to  avoid  some  problems  associated  with  the  spherical  harmonics.  Point 
masses  have  the  same  meaning  as  other  terms  used  in  the  geodetic  literature,  for 
example,  "buried  masses"  or  "mass  concentrations".  Point  masses  are  useful 
because  they  can  give  rise  to  the  disturbing  potential  from  which  geophysically 

(Received  for  publication  22  May  1979) 
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meaningful  quantities  such  as  geoid  undulations,  gravity  anomalies,  and  so  on, 
can  be  derived. 

The  disturbing  potential  added  to  the  normal  potential  results  in  the  total 

potential.  A  method  whereby  point  masses  give  rise  to  the  disturbing  potential  has 

been  presented  in  Reilly  and  Herbrechtsmeier.  1  Another  possibility  is  separation 

of  the  disturbing  potential  (T)  into  two  parts.  Tc  „  and  T„  „  ;  u  is  repre- 

n»  r •  ivl • 

sented  by  a  spherical  harmonic  expansion  and  Tp  M  by  the  point  mass  model. 

This  is,  in  fact,  the  avenue  pursued  in  this  report. 

The  idea  of  supplementing  the  spherical  harmonic  expansion  of  the  earth's 
potential  by  point  masses  is  relatively  recent.  The  purpose  for  introducing  the 
point  mass  parameters  into  an  adjustment  of  gravity  anomalies,  satellite  altimetry, 
and  perhaps  other  quantities  as  well,  is  to  add  fine  structure  to  a  geopotential 
model  based  on  spherical  harmonic  coefficients.  The  point  mass  locations,  includ¬ 
ing  their  depth,  theoretically  stipulated  with  respect  to  the  geoid  or  a  reference 
ellipsoid  are  the  object  of  judicial  choice,  while  the  point  mass  magnitudes  are  sub¬ 
jected  to  an  adjustment  in  a  generally  over-determined  system. 

A  first  adjustment  in  terms  of  spherical  harmonic  potential  coefficients  accom¬ 
modates  the  available  geoid  undulations,  obtained  from  satellite  altimetry,  as  well 
as  mean  gravity  anomaly  data.  A  subsequent  localized  adjustment  is  based  on  the 
altimetry  residuals  computed  in  a  region  of  interest  from  the  first  adjustment;  these 
residuals  are  accommodated  in  an  additional  data  fit  through  the  use  of  a  suitable 
number  of  new  parameters— the  point  mass  magnitudes.  The  point  mass  adjustment 
allows  a  relatively  small  number  of  parameters  to  express  local  details  of  the 
gravity  field  which  would  require  perhaps  a  hundredfold  increase  in  the  number  of 
parameters  if  the  adjustment  was  performed  only  on  the  basis  of  spherical  harmonics. 

The  results  are  contour  maps  of  geoid  undulations  and/or  gravity  anomalies 
constructed  from  values  computed  in  a  grid  extending  over  the  region  of  interest. 

No  point  masses  are  present  outside  the  local  region  so  that  the  geoid  there  essen¬ 
tially  coincides  with  the  global  geoid.  Contour  maps  of  standard  deviations  (sigmas) 
of  the  adjusted  and/or  predicted  quantities  can  also  be  constructed. 


2.  MATHEMATICAL  BACKGROUND 

In  the  spherical  approximation,  the  boundary  condition  and  Brun's  formula, 
respectively,  read 


1.  Reilly,  J.P.,  and  Herbrechtsmeier,  E.  H.  ( 1978)  A  systematic  approach  to 
modeling  the  geopotential  with  point  mass  anomalies,  J.  of  Geophys.  Res. 
83(No.  B2).  - ^ - 


A  g  =  (8T/8r) 


-  (2/R)  T. 


(1) 


r  =  R 

N  =  T/G  .  (2) 

where 

Ag  =  gravity  anomaly, 

N  =  geoid  undulation, 

T  =  disturbing  potential, 

R  =  earth's  mean  radius  (6371  km), 

G  =  average  value  of  gravity  (980  gals). 

These  two  formulas  can  be  found,  for  example,  on  pages  88  and  94  of  Heiskanen 
and  Moritz.  ^ 

Although  the  above  two  formulas  usually  refer  to  the  normal,  or  ellipsoidal, 
field,  in  this  report  they  indicate  a  higher  order  field,  specifically  the  spherical 
harmonic  (S.  H. )  model  after  the  first  adjustment.  Thus,  T  will  depict  the  disturb¬ 
ing  potential  associated  with  point  masses  (P.  M. )  and  will  describe  the  field  be¬ 
yond  the  S.  H.  model  (for  example,  beyond  the  14,  14  model).  Accordingly,  the 
quantities  A g,  N,  andT  appearing  in  Eqs.  (1)  and  (2)  could  be  written  with  the  sub¬ 
script  "tot"  if  they  refer  to  the  normal  field,  and  with  the  subscript  "P.  M.  "  if  they 
refer  to  a  higher  order  field  represented  by  the  S.  H.  model.  Since  it  holds  by 
definition  that 

w=U+Ttot. 

where 


W  =  actual  potential, 

U  =  normal  potential, 

Ttot  =  (total)  disturbing  potential, 

and  since  we  have  chosen  to  separate  the  total  disturbing  potential  into  two  parts. 


Ttot  =  TS.  H.  +  TP.  M.  * 


(3) 


2.  Heiskanen,  W.A. ,  and  Moritz,  H.  (1967)  Physical  Geodesy,  W.  H.  Freeman 
and  Co.,  San  Francisco,  CA. 
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Similar  relations  will  be  used  in  actual  adjustments,  not  only  for  predictions  but 
also  for  their  variances. 

The  parts  denoted  "S.  H.  "  can  be  treated  either  in  the  spherical  approximation 
or  in  a  more  accurate  form.  On  the  other  hand,  the  "P.  M. "  parts  can  always  be 
treated  in  the  spherical  approximation;  after  the  first,  S.  H.  adjustment,  the  resi¬ 
duals  to  be  accommodated  by  the  subsequent  P.  M.  adjustment  are  usually  at  least 
one  or  two  orders  of  magnitude  smaller  than  the  original  observations  (for  example, 
geoid  undulations)  so  that  this  kind  of  approximation  could  not  appreciably  contami¬ 
nate  the  resulting  geoid  representation,  gravity  anomalies,  and  so  on.  By  contrast, 
if  the  "total"  quantities  in  Eqs.  (4)  and  (5)  were  represented  by  a  model  containing 
the  spherical  approximation,  the  error  thus  introduced  could  amount  to  several 
decimeters  in  geoidal  heights.  It  is  thus  clear  that  if  the  P.  M.  model  containing 
the  spherical  approximation  were  used  in  an  adjustment  of  the  "total"  quantities 
referring  directly  to  the  normal  field,  a  decimeter  accuracy  in  the  geoid  determina¬ 
tion  could  not  even  be  approached,  due  to  this  approximation  alone. 

Since  we  shall  henceforth  deal  only  with  the  "P.  M.  "  part  of  the  quantities  in 
Eqs.  (3)  through  (5),  and  so  on,  the  subscript  "P.  M.  "  can  be  safely  omitted.  Based 
on  the  relation 

T.  =  £  (1/1..)  (kM>.  .  <6) 

i  j  ij  J 

where 

j.  =  distance  between  the  i-th  (observation)  point  and  the 

1J  j-th  point  mass, 

(kM)j  =  j-th  scaled  point  mass, 
from  Eqs.  (1)  and  (2)  we  can  derive 

Ag.  =  (1/R>£  (!/!..)  [<R2  -  Fy>/lJ  -  21  <kM).  .  (7) 
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(8) 


N  =  (l/G)  £  (III.  )  (kM).  . 

j  1J  J 

where 

l..  =  (R2  +  R2  -  2F..)1/2  , 

Rj  =  R  -  d  . 

F. .  =  RR,  cos  \l/. .  , 

1J  1 

cos  1^.^.  =  sin  sin  0^  +  cos  0.  cos  0^  cos  (Xj  -  X^.)  . 

The  spherical  approximation  becomes  further  apparent  from  the  fact  that  the  point 
masses  are  considered  to  be  located  on  a  sphere  of  radius  Rj,  at  the  depth  d  be¬ 
neath  the  surface  of  the  sphere  approximating  the  earth. 

3 

These  equations  are  adapted  from  formulas  which  appeared  in  Needham  and 
4  i 

Blaha  .  They  also  agree  with  the  algorithm  presented  in  Reilly  .  In  the  latter 

reference,  however,  the  formulas  referred  to  the  ellipsoid  and  not  to  a  higher  order 

surface;  with  the  exception  of  normal  gravity,  they  were  based  entirely  on  the 

spherical  approximation. 

3.  RESULTS  OF  ALTIMETER  DATA  REDUCTIONS 

The  data  in  the  first  global  adjustment  performed  at  AFGL  consist  of  the  avail¬ 
able  GEOS-3  satellite  altimeter  observations  and  the  global  set  of  1654  equal  5° 
mean  gravity  anomalies  supplied  by  R.  Rapp,  Ohio  State  University.  The  gravity 
anomaly  data  is  based  upon  approximately  38,000  1°  mean  anomalies.  Accuracy 
estimates  for  the  5°  mean  anomalies  are  depicted  in  Figure  1.  Of  the  1654  5°  mean 
anomalies  1507  were  based  directly  on  the  1°  anomalies,  while  the  remaining  147 
5°  means  were  obtained  by  interpolation.  The  mean  gravity  anomalies  serve  mainly 
to  make  the  spherical  harmonics  adjustment  possible.  Because  of  data  gaps  over 
the  continental  regions,  the  adjustment  could  not  be  performed  with  the  use  of 
satellite  altimetry  alone. 


3.  Needham,  P.  E.  (1970)  The  Formation  and  Evaluation  of  Detailed  Geopotential 

Models  Based  on  Point  Masses,  Department  of  Geodetic  Science. 

Report  No.  i49,  The  Ohio  State  University,  Columbus,  Ohio. 

4.  Blaha,  G.  (1977)  Refinements  in  the  Combined  Adjustment  of  Satellite  Altimetry 

and  Gravity  Anomaly  Data.  Report  of  DBA  Systems.  Inc.  :  AFGL-TK-77- 
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Figure  2  depicts  the  ground  tracks  of  1253  GEOS-3  passes,  over  oceanic 
regions,  used  in  this  global  adjustment  of  satellite  altimetry  and  gravity  anomalies. 
The  spherical  harmonic  model  is  complete  through  the  degree  and  order  (14,  14). 
The  altimeter  data  collected  over  the  North  Atlantic  region  are  further  used  in  a 
subsequent  point  mass  adjustment. 

In  the  second  or  residual  adjustment,  the  input  provided  by  satellite  altimetry 
completely  overrides  any  influence  by  the  gravity  anomaly  input  because  over  the 
region  of  interest,  such  as  the  North  Atlantic,  the  altimeter  data  density  is  hundreds 
or  thousands  of  times  higher  than  the  density  of  data  coming  from  the  anomalies 
and  the  weight  associated  with  each  altimeter  observation  is  considerably  higher 
than  the  weight  associated  with  any  given  gravity  anomaly.  Accordingly,  the 
parameters  in  the  second  adjustment  providing  the  local  geoidal  detail  are  the  point 
mass  magnitudes,  and  the  "observations"  in  the  area  of  concern  are  the  minus 
residuals  of  geoid  undulations  as  obtained  from  satellite  altimetry  in  the  first 
adjustment;  the  input  weights  are  unchanged  from  those  in  the  first  adjustment. 

The  results  from  the  first  adjustment  are  shown  in  Figures  3  and  4.  Figure  3 
depicts  contours  of  geoid  undulations  in  tO-m  intervals  and  Figure  4  depicts  con¬ 
tours  of  gravity  anomalies  in  10-mgal  intervals,  both  on  the  global  scale  over  the 
oceanic  regions.  The  reference  ellipsoid  used  in  these  computations  has  a  flatten¬ 
ing  of  1/298.  24  with  a  semi-major  axis  of  6,  378,  160  m  (Geodetic  Reference  System, 
1967). 

The  point  mass  locations  are  represented  by  dots  in  Figure  5.  Their  horizontal 
separation  is  2°  in  latitude  and  longitude,  forming  almost  equilateral  blocks  near 
the  equator.  A  2°  separation  is  suitable  for  representing  geoidal  detail  correspond¬ 
ing  to  a  half  wavelength  of  about  4°,  that  is,  two  parameters  per  half  wavelength  in 
each  direction.  The  depth  of  the  point  masses  corresponds  to  approximately  1.  6 
times  their  horizontal  separation,  in  this  case  about  350,  000  m.  This  configura¬ 
tion  has  been  confirmed  through  computer  simulations,  and  seems  to  be  a  suitable 
compromise  resulting  in  a  sufficiently  accurate  description  of  the  gravity  field, 
without  exceedingly  large  correlations  among  the  adjusted  parameters. 

Figure  5  also  shows  contours  of  geoid  undulations  in  5-m  intervals  obtained 
when  the  point  mass  adjustment  is  superimposed,  in  the  area  of  interest,  on  the 
previous  spherical  harmonic  adjustment.  One  can  detect  higher  geoidal  resolution 
within  the  area  covered  by  the  point  masses.  Remarkable  similarities  in  the  shape 
of  the  local  geoid  have  been  observed  between  this  geoidal  map  and  that  produced  in 
the  North  Atlantic  calibration  area  by  Marsh  and  Chang  in  1976,  who  based  their 
geoid  on  a  combination  of  5  ft  X  5  ft,  15  ft  X  15  ft,  and  1°  X  1°  surface  gravity  data 
and  the  GSFC  GM-8  Earth  Model. 
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Figure  5.  Geoid  Undulation  Contours  From  a  Spherical  Harmonic  14-Degree  Solution  and 
Point  Mass  Adjustment  Superimposed  on  the  Global  Adjustment  Using  Combinations  of 
Gravity  Anomaly  and  Altimetric  Data 


Finally,  Figure  6  depicts  the  predictions  of  free  air  gravity  from  the  same 
adjustment.  The  extreme  left-  and  right-hand  portion  of  this  figure  shows  the 
smooth  contour  lines  based  essentially  on  the  global  adjustment  in  terms  of  a 
(14,  14)  spherical  harmonic  model,  while  the  remaining  portion  shows  a  much  more 
detailed  (point  mass)  representation  of  the  gravity  field.  Such  a  detailed  repre¬ 
sentation  would  approximately  correspond  to  an  adjustment  in  terms  of  a  (45,  45) 
spherical  harmonic  model.  Irregular  behavior  of  a  point  mass  representation  has 
been  noticed  in  areas  containing  several  point  masses  but  very  few  observations, 
especially  when  depicting  gravity  anomalies,  which  represent  a  less  smooth  func¬ 
tion  than  the  geoid  undulations. 


4.  CONCLUSION 

In  many  present-day  adjustment  tasks  the  data  consist  of  mean  gravity  anomalies 
distributed  over  the  whole  globe  (with  some  exceptions)  and  of  densely  distributed 
satellite  altimeter  measurements  over  a  limited  area.  The  residuals  of  the  altim¬ 
eter  observations  obtained  in  a  global  adjustment  of  spherical  harmonic  potential 
coefficients  enter  the  subsequent,  localized  point  mass  adjustment. 

Good  results  in  terms  of  geoid  undulations  (both  adjusted  and  predicted)  and 
predicted  gravity  anomalies  can  be  obtained  if  the  shortest  half  wavelength  to  be 
represented  by  the  second  adjustment  corresponds  to  about  two  point  masses  in 
each  dimension  of  the  point  mass  grid.  The  depth  of  point  masses  which  has  been 
used  with  advantage  is  approximately  1.6  times  their  horizontal  separation. 

From  the  computer  simulations  documented  in  Blaha,  5  as  well  as  from  pre¬ 
liminary  data  reductions,  it  further  appears  that  in  order  to  arrive  at  reliable 
results,  the  point  mass  adjustment  should  be  based  on  a  fairly  high  degree  and 
order  spherical  harmonic  model  (for  example,  14,  14  model).  In  addition,  the 
geoidal  detail  to  be  represented  by  such  an  adjustment  should  be  contained  well 
within  the  point  mass  region.  This  region,  in  turn,  should  be  well  within  a  region 
of  high  density  observations. 


5.  Blaha,  G.  (1979)  Improved  Determination  of  the  Earth's  Gravity  Field,  Report 
of  Nova  University,  AFGL-Tr-79-0058. 
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ure  6.  Free  Air  Gravity  Anomaly  Contours  From  a  Spherical  Harmonic  14-Degree  Solution 
a  Point  Mass  Adjustment  Superimposed  on  the  Global  Adjustment  Using  Combinations  of 
vity  Anomaly  and  Altimetric  Data 


